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ABSTRACT 

The use of type Ic Super Luminous Supernovae (SLSN Ic) to examine the 
cosmological expansion introduces a new standard ruler with which to test theo¬ 
retical models. The sample suitable for this kind of work now includes 11 SLSNe 
Ic, which have thus far been used solely in tests involving ACDM. In this paper, 
we broaden the base of support for this new, important cosmic probe by using 
these observations to carry out a one-on-one comparison between the Rh = ct 
and ACDM cosmologies. We individually optimize the parameters in each cos¬ 
mological model by minimizing the statistic. We also carry out Monte Carlo 
simulations based on these current SLSN Ic measurements to estimate how large 
the sample would have to be in order to rule out either model at a ~ 99.7% con¬ 
fidence level. The currently available sample indicates a likelihood of ~70 — 80% 
that the Rh = ct Universe is the correct cosmology versus ~20 — 30% for the 
standard model. These results are suggestive, though not yet compelling, given 
the current limited number of SLSNe Ic. We find that if the real cosmology 
is ACDM, a sample of ~240 SLSNe Ic would be sufficient to rule out Rh = ct 
at this level of confidence, while ~480 SLSNe Ic would be required to rule out 
ACDM if the real Universe is instead Rh = ct. This difference in required sample 
size reflects the greater number of free parameters available to fit the data with 
ACDM. If such SLSNe Ic are commonly detected in the future, they could be a 
powerful tool for constraining the dark-energy equation of state in ACDM, and 
differentiating between this model and the Rh = ct Universe. 

Subject headings: cosmic background radiation - cosmological parameters - cosmology: 
observations - cosmology: theory - distance scale - supernovae: general 



1. Introduction 


A new type of luminous transient has been identified in recent years through the use 
of deep, wide surveys searching for supernovae in the local Universe (Quimby et al. 2005, 
2007; Quimby 2006; Smith et al. 2007; Drake et al. 2009; Rau et al. 2009; Kaiser et 
al. 2010; Gal-Yam 2012; Baltay et al. 2013). These super-luminous supernovae (SLSN) 
have peak magnitudes Mab < “21 mag and integrated burst energies 10^1 erg. They 
are therefore much brighter than both Type la SNe and the majority of core-collapse 
events. Their unusually high peak luminosities, hot blackbody temperatures and bright rest 
frame ultra-violet emission (which renders their continuum easily detectable at optical and 
near-infrared wavelengths at high redshifts; Cooke et al. 2012) allow them to be studied in 
concert with possible gamma-ray burst associations (Li et al. 2014; Cano Jakobsson 2014) 
and, more importantly, make them viable standardizable candles and distance indicators 
for use as cosmological probes (Inserra Sz Smartt 2014). 

Inserra et al. (2013) and Nicholl et al. (2013) used the classification term SLSN 
Ic to refer to all the hydrogen poor SNLSe, though there appear to be at least two 
observational groups in this category. These are distinguished via the terms SN2005ap-like 
and SN2007bi-like events, since these are the prototypes with the faster and slower evolving 
lightcurves. SLSNe Ic have now been discovered out to redshifts 2 : ~4 (Chomiuk et al. 
2011; Berger et al. 2012; Cooke et al. 2012; Howell et al. 2013), and appear to be rather 
homogeneous in their spectroscopic and photometric properties. 

The brighter events decline more slowly, not unlike the Phillips relation for Type 
la SNe, which raises the possibility of correlating their peak magnitudes to their decline 
over a fixed number of days to reduce the scatter (Rust 1974; Pskovskii 1977; Phillips 
1993; Hamuy et al. 1996). Without the use of such a relation, the uncorrected raw mean 
magnitudes show a dispersion of ~0.4 (e.g., Inserra & Smartt 2014). Correlating the peak 
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magnitude to the decline over 30 days reduces the scatter in standardized peak magnitudes 
to ±0.22 mag. And apparently using a magnitude-color evolution reduces this scatter even 
more, to a low value somewhere between ±0.08 mag and ±0.13 mag. 

It is therefore quite evident that SLSNe Ic may be useful cosmological probes, perhaps 
even out to redshifts much greater [z » 2) than those accessible using Type la SNe. 
The currently available sample, however, is still quite small; adequate data to extract 
correlations between empirical, observable quantities, such as lightcurve shape, color 
evolution and peak luminosity, are available only for tens of events. Our focus in this paper 
is specifically to study whether SLSNe Ic can be used—not only to optimize the parameters 
in ACDM (Li et al. 2014; Cano & Jakobsson 2014; Inserra & Smartt 2014), e.g., to refine 
the dark-energy equation of state but, also—to carry out comparative studies between 
competing cosmologies, such as ACDM versus the = ct Universe (Melia 2007; Melia & 
Abdelqader 2009; Melia & Shevchuk 2012; Melia 2013a). 

Like ACDM, the i7h = ct Universe is a Friedmann-Robertson-Walker cosmology that 
assumes the presence of dark energy, as well as matter and radiation. The principle 
difference between them is that the latter is also constrained by the equation of state 
p-f 3p = 0, in terms of the total pressure p and energy density p. In recent years, this model 
has generated some discussion concerning its fundamental basis, including claims that it 
is actually a vacuum solution, even though p ^ 0. However, all criticisms leveled against 
Rh = ct thus far appear to be based on incorrect assumptions and theoretical errors. A full 
accounting of this discussion may be found in Melia (2015) and references cited therein. 

In fact, the application of model selection tools in one-on-one comparisons between 
these two cosmologies has shown that the data tend to favor = ct over ACDM. Tests 
completed thus far include high-z quasars (Melia 2013b, 2014), gamma ray bursts (Wei et 
al. 2013), the use of cosmic chronometers (Melia Maier 2013) and, most recently, the 
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Type la SNe themselves (Wei et al. 2014). In all of these tests, information criteria show 
that Rh. = ct^ with the important additional constraint p + 3p = 0 on its equation of state, 
is favored over ACDM with a likelihood of ~ 90% versus only ~ 10%. 

Here, we broaden the comparison between i?h = ct and ACDM by now including 
SLSNe Ic in this study. In § 2 we briefly describe the currently available sample and our 
method of analysis. Our results are presented in § 3. We will find that the current catalog 
of SLSNe Ic suitable for this study already confirms the tendencies discussed above, though 
the statistics are not yet good enough to strongly differentiate between these two competing 
models. We show in § 4 how large the source catalog needs to be in order to rule out one or 
the other expansion scenario at a 3-sigma confidence level, and we present our conclusions 
in ? 5. 


2. Methodology 

Eleven of the SLSNe Ic identified by Inserra Sz Smartt (2014) are appropriate for this 
work, and we base our analysis on the methodology described in their paper and in Li et 
al. (2014). Briefly, the chosen SLSNe Ic must have well sampled light-curves around peak 
luminosity, and photometric coverage from several days pre-maximum to 30 days (in the 
rest frame) after the peak. This time delay appears to be optimal for use in the Phillips-like 
peak magnitude-decline relation (Inserra & Smartt 2014). 

All 11 of these events appear to be similar to the well-observed SN2010gx, and these 
decay rapidly after peak brightness. They belong to the group of 2005ap-like events, the 
first such SN discovered in this category. Other SLSNe Ic could not be included simply 
because of lack of sufficient temporal coverage, even though their identification has been 
securely classified in previous work (see, e.g., Leloudas et al. 2012; Quimby et al 2011; 
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Inserra et al. 2013; Chomiuk et al. 2011; Berger et al. 2012; Nicholl et al. 2014). 

The SNe listed in Table 1 were located in faint, dwarf galaxies, and were unlikely to 
have suffered significant extinction beyond the reddening induced by interstellar dust in 
our Galaxy. We here adopt the reddening corrections from Tables 1 and 2 of Inserra Sz 
Smartt (2014), who assumed a standard reddening curve with Ry = Ay/E{B — V) = 3.1. 
We also adopt their iiT-corrections and time dilation effects in order to obtain the absolute 
rest-frame peak magnitudes. This step is necessary due to the large (0.143 < 0 < 1.206) 
redshift coverage of the sample. These are listed in the table, along with the source names, 
their redshifts, apparent peak magnitudes (and filters), the magnitude decrease over 30 
days, and the color change between the 400nm and 520nm synthetic (restframe) bands at 
peak and 30 days later. 

In order to provide consistent, standarized comparative rest frame properties, the 
observed apparent magnitudes [m in Table 1) have been converted into defined, synthetic 
magnitudes. Since the SLSN Ic spectrum around 400nm is continuum dominated, Inserra 
Smartt (2014) defined a synthetic passband with an effective width of 80nm, centered 
at wavelength 400nm, having steep wings and a flat top. In Table 1, this is referred to as 
the 400nm band. All of the chosen events in this table have sufficient photometric coverage 
to allow the /^-correction to uniformly map the observed filter’s wavelength range to this 
400nm bandpass in the rest frame. The absolute magnitudes are then formally defined by 
the relation 

Afpeak(400) = m(400) - ^ 

= ruf- A'J!J^4oo - M , (1) 

where m(400) is the apparent magnitude mapped into the restframe 400nm band, ^ is the 
distance modulus calculated from the luminosity distance, rrif is the AB magnitude in the 
observed filter / {g, V, R, r, 0 , or z, as indicated in Table 1), Af is the Galactic extinction 



Table 1. Sample of SLSNe Ic 


SN 

z 

E{B-V) 

lUf 

Filter 


T^peak 

-^/^400 

m(400) 

AM3o(400) 

AM3o(400 -520) 

Refs. 

SN2011ke 

0.143 

0.01 

17.70 (g) 

9 - 

■> 400 

0.05 

-0.18 

17.83±0.20 

2.47±0.16 

0.54±0.17 

1, 8 

SN2012il 

0.175 

0.02 

18.00 (g) 

9 - 

■> 400 

0.09 

-0.18 

18.09±0.21 

2.19±0.16 

0.49±0.17 

1, 8 

PTFllrks 

0.190 

0.04 

19.13 (g) 

9 - 

■> 400 

0.17 

-0.19 

19.15±0.20 

2.62±0.14 

0.95±0.14 

1, 8 

SN2010gx 

0.230 

0.04 

18.43 (g) 

9 - 

■> 400 

0.13 

-0.23 

18.53±0.18 

2.00±0.19 

0.34±0.16 

2, 8 

SN2011kf 

0.245 

0.02 

18.60 (g) 

9 - 

■> 400 

0.09 

-0.13 

18.64±0.18 

1.49±0.16 

- 

1, 8 

LSQ12dlf 

0.255 

0.01 

18.78 (V) 

V - 

400 

0.03 

-0.27 

19.02±0.12 

l.OOiO.lO 

0.29±0.12 

3, 8 

PTF09cnd 

0.258 

0.03 

18.29 (R) 

R- 

400 

0.05 

-0.34 

18.58±0.22 

1.09±0.14 

- 

4, 8 

SN2013dg 

0.265 

0.01 

19.06 (g) 

9 - 

■> 400 

0.03 

-0.30 

19.33±0.20 

2.08±0.20 

0.51±0.17 

3, 8 

PSl-lObzj 

0.650 

0.01 

21.23 (r) 

r - 

^ 400 

0.02 

-0.58 

21.79±0.20 

1.70±0.14 

0.60±0.15 

5, 8 

PSl-lOky 

0.956 

0.03 

21.15 (i) 

i — 

>400 

0.06 

-0.73 

21.82±0.18 

1.31±0.15 

0.30±0.18 

6, 8 

SCP-06F6 

1.189 

0.01 

21.04 (z) 

z - 

> 400 

0.01 

-1.35 

22.38±0.20 

0.89±0.15 

- 

7, 8 


Notes: All values are from Inserra Sc Smartt (2014), except for m(400), which is calculated here. The error bars are directly 
from Figures 5 and 6 in Inserra Sc Smartt (2014). Columns; SN name; measured redshift; extinction; observed peak 
magnitude (AB system) and filter; Galactic extinction in the observed filter; calculation of the synthetic 400nm magnitude 
from the observed filter; apparent magnitude mapped into the 400nm passband, and its dispersion; magnitude decrease in 30 
restframe days and its dispersion; the color change between the 400nm and 520nm synthetic bands at peak and 30 days later, 
and its dispersion. References: (1) Inserra et al. (2013); (2) Pastorello et al. (2010); (3) Nicholl et al. (2014); (4) Quimby et 
al. (2011); (5) Lunnan et al. (2013); (6) Chomiuk et al. (2011); (7) Barbary et al. (2009); (8) Inserra Sc Smartt (2014). 







in the observed filter, and is the iiT-correction from the observed filter in Table 1 to 

the synthetic 400nm bandpass. The values of and m(400) are listed in Table 1. Note 
that m(400) is a cosmology-independent apparent magnitude. 

The peak magnitude-decline relation for SLSN Ic in the rest frame 400nm band is 
(Inserra & Smartt 2014) 

Mpeak(400) = Mo + aAM3o(400) , (2) 

where a is the slope, AM3o(400) is the decline at 30 days, and Mq is a constant representing 
the absolute peak magnitude at AM3o(400) = 0. Inserra & Smartt (2014) also found that 
Mpeak(400) appears to have a strong color dependence. Redder objects are fainter and also 
become redder faster. This peak magnitude-color evolution relation is given by 

Mpeak(400) = Mo + aAM3o(400 - 520) , (3) 

where AM3o(400 — 520) is the difference in color M(400) — M(520) between the peak and 
30 days later. (Note that Mo and a need not be the same in these two expressions.) 

With Equations (2) and (3), an effective (standardized) apparent magnitude may 
be obtained as follows: 

= m(400) — aAM^o . (4) 

The term aAMoo represents an adjustment due to the peak magnitude-decline relation, 
in terms of aAM3o(400), or the peak magnitude-color evolution relation, in terms of 
aAM3o(400 — 520), as the case may be. The effective apparent magnitude may also be 
expressed as (Perlmutter et ah 1997,1999) 

= T + 5 log^oWo ddz)] , (5) 

where Hq is the Hubble constant in units of km s~^ Mpc“^ and is the luminosity distance 
in units of Mpc. Here T is the “Mo-free” 400nm absolute peak magnitude, represented in 
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terms of the standardizable absolute magnitude Mq, according to the definition 

T = Mo-51ogio(i^o)+25 (6) 

(see Li et al. 2014 and references cited therein). 

The apparent correlations suggest the use of two “nuisance” parameters (a and T) 
whose optimization along with the model parameters decreases the overall scatter in the 
distance modulus. For each model, we therefore find the best fit by minimizing the 
statistic, defined as follows: 


X 


2 


E 


m,(400) - aAMso^i - T - 5logio[HodL{zi)]f 

<( 400 ) + 


( 7 ) 


In ACDM, the luminosity distance depends on several parameters, including 

Hq and the mass fractions = Pm/Pc? = Pr/Po a-nd fide = Pde/Pc? defined in terms of 
the current matter (pni)^ radiation (pr), and dark-energy (pde) densities, and the critical 
density pc = Sc^Hq/SttG. Assuming zero spatial curvature, so that flm + fir + fide = L the 
luminosity distance to redshift 2 : is given by the expression 

/* Z 

^ACDM(^) ^ C -\-z) [flm(l + zf + flr(l + z)^ + flde(l + z)3(l+^de)j ^ (g) 

-^0 Jo 

where pde = '^dePde is the dark-energy equation of state. The Hubble constant Hq cancels 
out in Equation (7) when we multiply di by Hq, so the essential remaining parameters 
in flat ACDM are fl^ and u^de- If we further assume that dark energy is a cosmological 
constant with Wde = “I 5 then only the parameter fl^ is available to fit the data. 


In the = ct Universe (Melia 2007; Melia Sz Abdelqader 2009; Melia & Shevchuk 
2012), the luminosity distance depends only on Hq, but since here too the Hubble constant 
cancels out in the product Hddi, there are actually no free (model) parameters left to fit 
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the SLSN Ic data. In this cosmology, 

d2'^=^\z) = ^il + z)Hl + z). (9) 

-no 

3. Results 

We have assumed that SLSNe Ic can be used as standardizable candles and applied the 
AM 30 decline relation (with 11 objects) and the peak magnitude-color evolution relation 
(with 8 objects) to compare the standard (ACDM) model with the Rh = ct Universe. In 
this section, we discuss how the fits have been optimized, first for ACDM, and then for 
Rh = ct. The outcome for each model is more fully described and discussed in subsequent 
sections. 


3.1. ACDM 

In the most basic ACDM model, the dark-energy equation of state parameter, 
is exactly —1. The Hubble constant Hq cancels out in Equation (7) when we multiply 
by Hq, so the essential remaining parameter is Type la SN measurements (see, e.g., 
Perlmutter et al. 1998, 1999; Riess et al. 1998; Schmidt et al. 1998), CMB anisotropy 
data (e.g., Hinshaw et al. 2013), and baryon acoustic oscillation (BAO) peak length scale 
estimates (e.g., Samushia Ratra 2009), strongly suggest that we live in a spatially flat, 
dark energy-dominated universe with concordance parameter values ~ 0.27 and Hq ^ 70 
km s“^ Mpc“^. 

We will therefore first attempt to fit the data with this concordance model, using 
prior values for all the model parameters, but not the two “nuisance” parameters a and T. 
For the AM 30 decline relation (with 11 objects), the resulting constraints on a and T are 
shown in Figure 1(a). For this fit, we obtain a = 0 . 57 ^ 0^7 ^ = “6.75^032 
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Fig. 1.— l-3a constraints on a and T for the concordance model, using the AAfso decline 
relation (panel a) and the peak magnitude-color evolution relation (panel b). 

a degree of freedom of Xdof ~ 27.04/9 = 3.00, remembering that all of the ACDM 

parameters are assumed to have prior values, except for the “nuisance” parameters a and 
T. (The corresponding data and best fit are shown in the top lefthand panel of Figure 4.) 
For the peak magnitude-color evolution relation (with 8 objects), the resulting constraints 
on a and T are shown in Figure 1(b). The best-fit parameters are a = 2 . 00 ^Q 7 g (Icr), 

T = -6.65ig;^? (Ifj). The degree of freedom for the concordance model with these 

optimized “nuisance” parameters is Xdof “ 2.48/6 = 0.41 (see also the top righthand panel 
in Figure 4). 

If we relax some of the priors, and allow fo be a free parameter, we obtain the 1 - 3(7 
constraint contours shown in the top panel of Figure 2 for the AA /30 decline relation, and 
the bottom panel for the peak magnitude-color evolution relation. These contours show 
that at the la level, the optimized parameter values are a = 0.46^Q}g, T = —6.41 /q 33, and 

= 0 . 62 / 334 . We find that the x^ P^^ degree of freedom for the optimized flat ACDM 
model is Xdof ~ 24.89/8 = 3.11. The bottom panel of Figure 2 shows the 1 - 3(7 constraints 
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-5.5 






Fig. 2. — Top panel: Constraints on a, T and Hni for the flat ACDM model, using the AM 30 
decline relation. Bottom panel: Same as the top panel, but using the peak magnitude-color 
evolution relation. 
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on a, T, and Qrn for the flat ACDM model, using the peak magnitude-color evolution 
relation. The contours show that at the Icr-level, a = 2.05^094, T = —6.52lQQg, but that 
r^ni is poorly constrained; only a lower limit of ~ 0.05 can be set at this confidence level 
and the best-fit is 0.59. The por degree of freedom is Xdof ~ 2.15/5 = 0.43. 


3.2. The = ct Universe 

The Rh. = ct Universe has only one free parameter, Hq, but since the Hubble constant 
cancels out in the product there are actually no free (model) parameters left to fit 

the SLSN Ic data. The results of fitting the AM 30 decline relation with this cosmology 
are shown in Figure 3(a). We see here that the best fit corresponds to a = O.bO/g'^^ 

(Icr) and T = —O.dSto^si (1*^)- With 11 — 2 = 9 degrees of freedom, the reduced is 
X^^^ = 25.79/9 = 2.87. The results of fitting the peak magnitude-color evolution relation 
with this cosmology are shown in Figure 3(b). We see here that the best fit corresponds 
to a = 1 . 88/973 (f^) ^ ~ (1(^)- With 8 — 2 = 6 degrees of freedom, the 

reduced x^ is Xdoi = 2.19/6 = 0.37. 

To facilitate a direct comparison between ACDM and i?h = ct^ we show in Figure 4 the 
Hubble diagrams for SLSNe Ic. In the left panel of Figure 4, the effective magnitudes 
of 11 SLSN Ic are plotted as solid points, together with the best-fit theoretical curves (from 
top to bottom) for the concordance model (with prior values of the parameters, and with 
a = 0.57 and T = —6.75), for the optimized flat ACDM model (with Dm = 0.62, and with 
a = 0.46, and T = —6.41). and for the = ct Universe (with a = 0.50 and T = —6.48). 
The Hubble diagrams derived using the peak magnitude-color evolution relation are shown 
in the righthand panels of Figure 4. Here, the effective magnitudes of 8 SLSN Ic are 
plotted as solid points, together with the best-fit theoretical curves for the concordance 
model (with prior values of the parameters, and with a = 2.00 and T = —6.65), for the 
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Fig. 3.— l-3cr constraints on a and T for the = ct Universe, using the AM 30 decline 
relation (panel a) and the peak magnitude-color evolution relation (panel b). 

optimized flat ACDM model (with Qm = 0.59, and with a = 2.05, and T = —6.52), and for 
the Rh = ct Universe (with a = 1.88 and T = —6.44). Strictly based on their Xdof values, 
the optimized ACDM model and the i?h = ct Universe appear to fit the SLSN Ic data 
comparably well. However, because these models formulate their observables (such as the 
luminosity distances in Equations 8 and 9) differently, and because they do not have the 
same number of free parameters, a comparison of the likelihoods for either being closer to 
the ‘true’ model must be based on model selection tools. 


3.3. Model Selection Tools 

Several information criteria commonly used in cosmology (see, e.g., Melia Sz Maier 2013, 
and references cited therein) include the Akaike Information Criterion, AIC = x^ + 2n, where 
n is the number of free parameters (Akaike 1973; Takeuchi 2000; Liddle 2004, 2007; Tan & 
Biswas 2012), the Kullback Information Criterion, KIC = + 3n (Bhansali & Downnham 

1977; Cavanaugh 1999, 2004), and the Bayes Information Criterion, BIC = + {lnN)n^ 
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Fig. 4.— Hubble diagrams for SLSN Ic. The effective apparent magnitudes of SLSN Ic are 
plotted as solid points, together with their corresponding best-fit theoretical curves. Left: 
The effective magnitudes are calculated using the AM^q decline relation. Right: The effective 
magnitudes are calculated using the peak magnitude-color evolution relation. 
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where N is the number of data points (Schwarz 1978; Liddle et ah 2006; Liddle 2007). 
With AICq characterizing model A4a, the unnormalized confidence that this model is true 
is the Akaike weight exp(—AICa/2). Model AIq has likelihood 


P(Ma) 


exp(-AICa/2) 

exp(-AICi/2) +exp(-AIC 2 / 2 ) 


( 10 ) 


of being the correct choice in this one-on-one comparison. The difference AAIC = 
AIC 2 — AlCi determines the extent to which Mi is favoured over A^ 2 - For Kullback and 
Bayes, the likelihoods are defined analogously. 


With the optimized fits we have reported in this paper, our analysis of the AM 30 
decline relation (with 11 objects) shows that i?h = ct is favoured over the flat ACDM model 
with a likelihood of 5 =^^ 63.4% versus 36.6% using AIC, 74.1% versus 25.9% using KIC, 
and 67.9% versus 32.1% using BIC. In our one-on-one comparison using the peak 
magnitude-color evolution relation (with 8 objects), the = ct Universe is preferred over 
ACDM with a likelihood of ^ 72.7% versus 27.3% using AIC, ~ 81.5% versus ~ 18.5% 
using KIC, and 73.5% versus 26.5% using BIC. 


4. Monte Carlo Simulations with a Mock Sample 

These results are interesting, though the current sample of SLSNe Ic is clearly too 
small for either model to be ruled out just yet. However, this situation will change with the 
discovery of new SNe, particularly at redshifts 2 > 2. To anticipate how well the SLSN Ic 
catalog obeying the peak magnitude-color evolution relation may be used to constrain the 
dark-energy equation of state in ACDM, and to differentiate between this standard model 
and the 7?h = ct Universe, we will here produce mock samples of SLSNe Ic based on the 
current measurement accuracy. In using the model selection tools, the outcome A = AICi — 
AIC 2 (and analogously for KIC and BIC) is judged ‘positive’ in the range A = 2 — 6 , 
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‘strong’ for A = 6 — 10, and ‘very strong’ for A > 10. In this section, we will estimate the 
sample size required to significantly strengthen the evidence in favour of = ct or ACDM, 
by conservatively seeking an outcome even beyond A ~ 11.62, i.e., we will see what is 
required to produce a likelihood ^ 99.7% versus ^ 0.3%, corresponding to a 3cr confidence 
level. 

We will consider two cases: one in which the background cosmology is assumed to be 
ACDM, and a second in which it is i?h = ct, and we will attempt to estimate the number 
of SLSNe Ic required in each case to rule out the alternative (presumably incorrect) model 
at a ~ 99.7% confidence level. The synthetic SLSNe Ic are each characterized by a set of 
parameters denoted as ( 2 , m[400], AM3o[400 — 520]). We generate the synthetic sample 
using the following procedure: 

1 . Since a subclass of broad-lined Type Ic SNe are observed to be associated with 
Gamma-ray bursts (GRBs; Hjorth et al. 2003; Stanek et al. 2003), we simulate the 
2 -distribution of our sample based on the observed 2 -distribution of GRBs. Shao et al. 
(2011) found that the 2 -distribution of GRBs appears to be asymptotic to a parameter-free 
probability density distribution, f{z) = ze~^. The redshift 2 of our SLSN Ic events is 
generated randomly from this function. Since they can be discovered out to redshifts 2^4 
(e.g., Chomiuk et al. 2011; Berger et al. 2012; Cooke et al. 2012; Howell et al. 2013), the 
range of redshifts for our analysis is [0,4]. We assign the absolute peak magnitude M(400) 
uniformly between -23.0 and -21.0 mag, based on the current SLSN Ic measurements. 

2. The synthetic apparent peak magnitude m(400) is calculated using the relation 
m(400) = M(400) + fJ-{z), where /i( 2 ) is the distance modulus at 2 , for a flat ACDM 
cosmology with = 0.27 and Hq = 70 km s“^ Mpc“^ (§ 4.1), or the Rh = ct Universe 
with Hq = 70 km s“^ Mpc“^ (§ 4.2). 


3. With the selected 2 and M(400) values, we first infer a AM3o(400 — 520) from the 
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peak magnitude-color evolution relation in a flat ACDM cosmology with = 0.27 and 
Hq = 70 km Mpc“^ (§ 4.1), or the _Rh = ct Universe with Hq = 70 km s“^ Mpc“^ (§ 4.2). 
And then we add scatter to this relation by assigning a dispersion to the AM3o(400 — 520) 
value; i.e., we randomly map this quantity to the new value AM3 q( 400 — 520) assuming a 
normal distribution with a center at AM3o(400 — 520) and a dispersion a = 0.1 mag. This 
value of a is typical for the current (observed) peak magnitude-color evolution relation, 
which yields a standard deviation a = 0.08 mag for the linear fit (Inserra Sz Smartt 2014). 
If AM3 q(400 — 520) > 0.0 mag, this SLSN Ic is included in our sample. Otherwise, it is 
excluded. 

4. We next assign “observational” errors to m(400) and AM3 q( 400 — 520). We will 
assign a dispersion ai = 0.05 + 0.05^^ to each event i, where ki is a random number between 
0 and 1. 

This sequence of steps is repeated for each new SLSN Ic in the sample, until we reach 
the likelihood criterion discussed above. As with the real 11-SLSN Ic sample, we optimize 
the model fits by minimizing the function in Equation (7). 

4.1. Assuming ACDM as the Background Cosmology 

We have found that a sample of at least 240 SLSNe Ic is required in order to rule out 
i?h = ct at the ~ 99.7% confidence level. The optimized parameters corresponding to the 
best-fit rcCDM model for these simulated data are displayed in Figure 5. To allow for the 
greatest flexibility in this fit, we relax the assumption that dark energy is a cosmological 
constant with Wde = “1, and allow Wde to be a free parameter, along with flm. Figure 5 
shows the 1-D probability distribution for each parameter (flmi T), and 2-D plots of 

the l-3a confidence regions for two-parameter combinations. The best-fit values for rcCDM 
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using the simulated sample with 240 SNe in the ACDM model are = 0.28loo4 
Wde = « = 2.2n°:}° (la), and T = -6.80i°:}° (la). 

To gauge the impact of these constraints more clearly, we show in Figure 6 the 
confidence regions (shaded, with red contours) for a^nd Wde using these 240 simulated 
SLSNe Ic (the same as the bottom left-hand panel of Figure 5), and compare these to the 
constraint contours for the 580 Union2.1 Type la SN data (Suzuki et al. 2012) (represented 
by the blue contours in Figure 6). It is straightforward to see how effectively the SLSNe Ic 
could be used as a cosmological tool, because the confidence regions resulting from their 
analysis are smaller (and narrower for flm) than those corresponding to the Type la events. 
The better constraints are mainly due to the fact that SLSNe Ic are distributed over a much 
wider redshift range, extending towards high-z, where tighter constraints on the model can 
be achieved. 

In Figure 7, we show the corresponding 2-D contours in the a — T plane for the 
7?h = ct Universe. The best-fit values for the simulated sample are a = 2.22lQQg (la) and 
T = -6.53l°:°l (la). 

Since the number N of data points in the sample is now much greater than one, the 
most appropriate information criterion to use is the BIG. The logarithmic penalty in this 
model selection tool strongly suppresses overfitting if N is large (the situation we have 
here, which is deep in the asymptotic regime). With N = 240, our analysis of the simulated 
sample shows that the BIG would favor the rcGDM model over = ct by an overwhelming 
likelihood of 99.7% versus only 0.3% (i.e., the prescribed 3a confidence limit). 
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Fig. 5. The 1-D probability distributions and 2-D regions with the l-3cr contours cor¬ 
responding to the parameters and T in the best-fit rcCDM model, using the 

simulated sample with 240 SLSNe Ic, assuming ACDM as the background cosmology. 
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4.2. Assuming i?h = ct as the Background Cosmology 

In this case, we assume that the background cosmology is the i?h = ct Universe, and 
seek the minimum sample size to rule out 'u;CDM at the 3(7 confidence level. We have 
found that a minimum of 480 SLSNe Ic are required to achieve this goal. To allow for the 
greatest flexibility in the rcCDM fit, here too we relax the assumption of dark energy as 
a cosmological constant with Wde = “U and allow to be a free parameter, along with 
In Figure 8, we show the 1-D probability distribution for each parameter 
a, T), and 2-D plots of the l-Scr confidence regions for two-parameter combinations. The 
best-fit values for rcCDM using this simulated sample with 480 SLSNe Ic are = 0.0, 
Wde = —0.33io;o6 (I*^)’ ^ = 2.05io;o7 (I<^)’ ^ ^ “6.52io;o5 (I<^)- Note that the 

simulated SLSNe Ic give a good constraint on Wde, but a weak constraint on only an 
upper limit of 0.09 can be set at the la confidence level. 

The corresponding 2-D contours in the a — T plane for the = ct Universe are 
shown in Figure 9. The best-fit values for the simulated sample are a = 2.05ioo6 (I^) 

T = “6*52^0 02 (I*^)- These are similar to those in the standard model, but not exactly 
the same, reaffirming the importance of reducing the data separately for each model being 
tested. With N = 480, our analysis of the simulated sample shows that in this case the BIG 
would favor R]^ = ct over rcCDM by an overwhelming likelihood of 99.7% versus only 0.3% 
(i.e., the prescribed 3a confidence limit). 


5. Conclusions 


It is quite evident that SLSNe Ic may be useful cosmological probes, perhaps even out 
to redshifts much greater {z >> 2) than those accessible using Type la SNe. The currently 
available sample, however, is still quite small; adequate data to extract correlations 
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between empirical, observable quantities, such as lightcurve shape, color evolution and peak 
luminosity, are available only for tens of events. In this paper, we have proposed to use 
SLSNe Ic for an actual one-on-one comparison between competing cosmological models. 
This must be done because the results we have presented here already indicate a strong 
likelihood of being able to discriminate between models such as ACDM and = ct. Such 
comparisons have already been made using, e.g., cosmic chronometers (Melia & Maier 
2013), Gamma-ray bursts (Wei et al. 2013), and Type la SNe (Wei et al. 2015). 

We have individually optimized the parameters in each model by minimizing the 
statistic. With the optimized fits we have reported above, our analysis of the AM 30 decline 
relation (with 11 objects) shows that i?h = ct is favoured over the flat ACDM model with a 
likelihood of ^ 63 — 74% versus 26 — 37% (depending on the information criterion). In our 
one-on-one comparison using the peak magnitude-color evolution relation (with 8 objects), 
the i?h = ct Universe is preferred over ACDM with a likelihood of 73 — 82% versus 
18 - 27%. 

But though SLSN Ic observations currently tend to favor = ct over ACDM, the 
known sample of such measurements is still too small for us to completely rule out either 
model. We have therefore considered two synthetic samples with characteristics similar to 
those of the 8 known SLSN Ic measurements, one based on a ACDM background cosmology, 
the other on R]^ = ct. From the analysis of these simulated SLSNe Ic, we have estimated 
that a sample of about 240 such events are necessary to rule out i?h = ct at a ~ 99.7% 
confidence level if the real cosmology is in fact ACDM, while a sample of at least 480 SNe 
would be needed to similarly rule out u^CDM if the background cosmology were instead 
= ct. The difference in required sample size results from rcCDM’s greater flexibility in 
fitting the data, since it has a larger number of free parameters. 

Our simulations have also shown that a moderate sample size of ^ 250 events could 
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reach much tighter constraints on the dark-energy equation of state Wde and on the matter 
density fraction Qrn than are currently available with the 580 Union2.1 Type la SNe. If 
SLSNe Ic can be commonly detected in the future, they have the potential of greatly 
refining the measurement of cosmological parameters, particularly the dark-energy equation 
of state Wde- 

Assembling samples of this size may be feasible with upcoming surveys. For example, 
the planned survey SUDSS (Survey Using Decam for Super-luminous Supernovae) using 
the Dark Energy Camera on the CTIO Blanco 4m telescope Inserra Sz Smartt (2014), and 
the Subaru/Hyper Suprime-Cam deep survey (Tanaka et al. 2012), have a goal of discovery 
several hundred SLSNe out to z ^ 4 over 3 years by imaging tens of square degrees to 
a limiting magnitude ~ 25 every 2 weeks. These numbers are based on an estimated 
production in the local universe (Quimby et al. 2013), who reported a rate of 32^26 events 
Gpc“^ yr“^ at a weighted redshift of (z) = 0.17. This number is low (0.01%) compared 
to that of Type la SNe. McCrum et al. (2015) estimate a Type Ic SN rate ~ 10“^ of the 
overall core-collapse SN rate within 0.3 < z < 1.4, though this number could be higher 
at 2 > 1.5 (Cooke et al. 2012) due, perhaps, to a decreasing metallicity. Moreover, the 
increase in cosmic star formation rate would boost the absolute numbers of SLSNe. As far 
as observations from the ground are concerned, assembling a sample of several hundred 
Type Ic SLSNe over 3 to 5 years therefore looks quite promising, particulary with the 
surveying capability of LSST (see, e.g., Lien & Fields 2009), which should cover the whole 
sky every 2 nights, down to a limiting magnitude ~ 24. The expected rate of discovery of 
core-collapse SNe with this survey is expected to be ^ 1 — 2 s“^ out to a redshift z ~ 2. 
Given that ~ 10“^ of these are expected to be Type Ic SNe in this redshift range, one 
should expect LSST to assemble a ~ 500 SLSN sample in less than a year. The situation 
from space is even more exciting. Detecting Type Ic SLSNe out to redshifts z ^ 10 via their 
restframe 400nm and 520nm bands is plausible with, e.g., EUCLID (Laureijs et al. 2011), 
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WFIRST and the James Webb Space Telescope (specifically the NIRCam)!^ 

A major problem with this approach right now, however, is that one must rely on the 
use of a luminosity-limited sample of supernovae, i.e. those selected to be super-luminous 
with peak magnitudes Mab < “21 mag and integrated burst energies erg (Inserra & 

Smartt 2014), in order to test the luminosity distances. The difficulty is that the luminosity- 
limited sample may turn out to be different with different background cosmologies, given 
that their inferred luminosities are themselves dependent on the models. This situation is 
unlikely to change as the sample grows, so it would be necessary to find a way of identifying 
these SNe, other than simply through their magnitudes. Fortunately, in the case of Rh = d 
versus ACDM, even though their luminosity distances are formulated differently, it turns 
out that the ensuing distance measures derived from these are quite similar all the way 
out to z ~ 6. As such, Type Ic SLSNe that make the cut for one model and not the other 
are the exception rather than the rule. In other words, the luminosity-limit applied to the 
sample examined here does not bias either model very much. But this may not be true in 
general, and an alternative method of selection is highly desirable. 
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Fig. 6.— One-cr, 2-a and 3-cr confidence regions for the icCDM model using simulated SLSNe 
Ic (shaded, with red contours), compared with those (dashed blue contours) associated with 
the 580 Union2.1 Type la SNe. 
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Fig. 7.— The 2-D region with the 1-3(7 contours for the parameters a and T in the i?h = 
ct Universe, using a sample of 240 SLSN Ic, simulated with ACDM as the background 
cosmology. The simulated model parameters were = 0.27, = 0.73, and Hq = 70 km 

s-^ Mpc-h 
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Fig. 8.— Same as Figure 5, except now with = ct as the (assumed) background cosmology. 
The simulated model parameter was Hq = 70 km s~^ Mpc“^. 
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Fig. 9.— Same as Figure 7, except now with i?h = ct as the (assumed) background cosmology. 
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